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Abstract. By a proper choice of the excitation energy per nucleon we analyze the 
mass distributions of the nuclear fragmentation at various excitation energies. Starting 
from low energies (between 0.1 and 1 MeV/nucleon) up to higher energies about 12 
MeV/n, we classified the mass yield characteristics for heavy nuclei {A > 200) on the 
basis of Statistical Multifragmentation Model. The evaluation of fragment distribution 
with the excitation energy show that the present results exhibit the same trend as the 
experimental ones. 
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1. Introduction 

Properties of nuclear reactions have been under investigation for several decays. 
Experimental and theoretical studies on nuclear reactions are very important not only 
for the context of nuclear physics but also for our understanding of astrophysical events, 
such as supernova explosion and formation of neutron stars |1]. At low excitation 
energies up to 1 MeV per nucleon, the fission of compound nucleus or its evaporation 
to small particles have been observed both theoretically and experimentally. At the 
excitation energies in between 1 and 3 MeV per nucleon one may observe a U-shape 
distribution corresponding to partitions with a few small fragments and one big residual 
fragment, which looks like an evaporative emission. At excitation energies greater than 
3 MeV per nucleon one may observe a mass distribution of intermediate mass fragments 
(nuclear multifragmentation). During this process, the hot and compressed nuclei tend 
to expand in thermodynamical equilibrium. Then they enter the region of subsaturation 
densities, where they become unstable to density fluctuations and break up into the 
fragments. In this case, it is believed that a liquid-gas phase transition is manifested 
[2], [31 m E] . At higher energies 8 — 9 MeV/n, big fragments disappear, and an exponential 
fall of the mass distribution with mass number A may be seen. 



Mass distributions for nuclear disintegration from fission to evaporation 



2 



By a proper choice of the excitation energy per nucleon we analyze the mass 
distributions of the nuclear fission of heavy elements such as ^ss-q g^j^^ 226j^g^ various 
excitation energies on the basis of Statistical Multifragmentation Model (SMM). We 
have also discussed the mass distribution of the nuclear fragmentation at higher energies 
up to 12 MeV/n. 

2. Calculations and results 

Nuclear fission is one of the most interesting channel of de-excitation of heavy nuclei at 
low excitation energies up to 1 MeV/nucleon, at which nuclear density exhibits small 
fluctuations around the equilibrium density 0.16 fm~^. In order to study compound 
nucleus fission one should consider a correct physical description in the light of liquid 
drop model with shell effects at various deformation modes, where the determination 
of angular momentum and dynamical variables becomes very important. In our 
calculations, we consider the Bohr-Wheeler approach [6] in SMM, where the partial 
width of the compound nucleus is given by 



where po is the level density of compound nucleus, pf the level density at the saddle 
point, E* the excitation energy for the nucleus, and Bf the height of the fission barrier. 
Here, the shell effect on the level densities of the nucleus can be neglected at high 
energies [7]. The fission barrier is determined by Myers and Swiatecki [8], and we use 
the results in Ref. [9] for the level density at saddle point. 

For orientation, we have shown the results of our calculations at low excitation 
energies for the fission of 238u in Fig. 1. In this figure, one may see a double peaked 
distribution at an excitation energy of 0.1 MeV/n on the upper-left panel, which can be 
interpreted as asymmetric fission, where the mass ratio of the most probable fragments 
is about 3 to 2. In the upper-right panel one may see a triple peaked mass distribution, 
which is usually interpreted in terms of a single symmetric peak and a double-peaked 
distributions. In the lower panels we show the characteristic of symmetric fission, 
associated with higher excitation energies. The evaluation of mass yield distribution 
with the excitation energy show that the present results are in agreement with those 
obtained experimentally in the reactions induced by protons and deuterons on ^^^Ra 



It is also observed that the proportion of the symmetric component increases with 
increasing excitation energy, in agreement with experimental results. In order to make 
a comparison with p!Q], we have chosen the nucleus ^^^Ra, and showed the results of 
our calculations in Fig. 2, from which one may see a triple peaked mass distribution 
at 0.1 MeV/n, and a single symmetric distribution at 0.2 MeV/n. The triple peaked 
distribution can be interpreted in terms of a symmetric fission associated with higher 
excitation energies, and an asymmetric component associated with lower excitation 
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Figure 1. Mass yield distributions for the fission of 238]j various excitation 
energies between 0.1 and 0.4 MeV per nucleon within SMM. 
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Figure 2. The left panel shows the mass yields of the fission of Ra as a triple 
peaked mass distribution at 0.1 MeV/n (denoted by the open circles), and a 
single symmetric distribution at 0.2 MeV/n (the full circles). The right panel 
illustrates the percentage of the same distributions. 

energies. In the right top and bottom panels, we also show the results in percentage of 
the variation of the mass yield with mass number A. These results are very similar to 
the obtained experimentally in Ref . [10] . We also review the fragmentation phenomena 
at higher excitation energies on the basis of SMM. According to SMM, one assumes a 
micro-canonical ensemble of breakup channels, and the system should obey the laws of 
conservation of energy E*, mass number A and charge number Z. The probability of 
generating any breakup channels is assumed to be proportional to its statistical weight 
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as 

Wj (xexp{Sj{E*,A,Z)) (2) 

where Sj denotes the entropy of a multifragment state of the breakup channel j. The 
breakup channels are generated by Monte Carlo method according to their weights. 
Light fragments with mass number and charge number are considered as stable 
particles (nuclear gas) with masses and spins taken from the nuclear tables. Only 
translational degrees of freedom of these particles contribute to the entropy of the 
system. Fragments are treated as heated nuclear liquid drops, and their individual 
free energies are parameterized as a sum of the bulk, surface. Coulomb and symmetry 
energy contributions 

Fa,z = F^,z + Fa,z + Ea,z + Ea,z- (3) 

The bulk contribution is given by F^^ — (^o — ^^/^o)^, where T is the temperature, 
the parameter Eq is related to the level density, and Wq = 16 MeV is the binding 
energy of infinite nuclear matter. Contribution of the surface energy is F^^ = 
BoA^/^KT^ - T^)/{T^ + T^)f/\ where = 18 McV is the surface coefficient,' and 
Tc = 18 MeV the critical temperature of the infinite nuclear matter. Coulomb energy 
contribution is = cZ'^/A^^^, where c denotes the Coulomb parameter obtained in 
the Wigner-Seitz approximation, c = (3/5)(e^/ro)(l — (p/po)^^"^), with the charge unit 
e. To = 1.17 fm, and po is the normal nuclear matter density (0.15 fm~^). And finally, 
the symmetry term is E^'^ — 'y{A — 2Zy/A, where 7 = 25 MeV is the symmetry 
energy parameter. All the parameters given above are taken from the Bethe-Weizscker 
formula and correspond to the assumption of isolated fragments with normal density in 
the freeze-out configuration. 
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Figure 3. Mass yield distributions for the fragmentation of ^^^Ra nucleus, at 
excitation energies of E* = 4, 5 and 8 MeV/n, within SMM. 

In Fig. 3 we show the results of our calculation for typical mass distributions for the 
fragmentation of ^^^Ra nucleus, at excitation energies of E* = 4, 5 and 8 MeV/n. One 
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may see from this figure that at an excitation energy of E* = 4 MeV/n (corresponding 
temperature T < 5 MeV), there is a U-shape distribution corresponding to partitions 
with few small fragments and one big residual fragment. In the so called transition region 
(T ?5i 5 — 6 MeV), however, variation of T with E* exhibits a plateau-like behavior that 
can be interpreted as a sign of liquid-gas phase transition in the system [HI [121 fT3] . 
and one observes a smooth transformation for E* = 5 MeV/n in the same figure. At 
high temperatures (T > 6 MeV) for E* = 8 MeV/n, the big fragments disappear and an 
exponential-like fall-off is observed. For higher excitation energies (E* > 12 MeV/n) the 
nuclei tend to evaporate to nucleons and small fragments. All these results are in good 
agreement with experimental data [5l [HI UHl |16l |17] . In view of these theoretical results 
it is instructive to demonstrate the possibility of the application of such approaches for 
the analysis of experimental data for nuclear reactions and astrophysical studies [Tl 118]. 
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